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Abstract 

Possible effects due to nonequilibrium dynamics in the Affleck-Dine mechanism of baryogene- 
sis are examined. Using the closed-time-path formalism, the quantum fluctuation and the back 
reaction of the Affleck-Dine scalar field are incorporated self-consistently into the dynamical equa- 
tions of the system by invoking a nonperturbative Hartree approximation. It is found that such 
nonequilibrium effects can significantly affect the amount of baryon asymmetry that can be gener- 
ated. In particular, it is possible to generate the observed baryon asymmetry with suitable initial 
conditions. The methodology described in this paper as well as some of the results obtained are 
quite general, and can be applied to any complex scalar field in a cosmological background. 
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I. INTRODUCTION 



It is evident that no concentration of antimatter exists within the Solar system and the 
Milky Way. The absence of annihilation radiation from the Virgo cluster indicates that 
antimatter can hardly be found within a 20 Mpc scale. A study of the contribution of 
annihilation radiation near the matter-antimatter boundaries to the cosmic diffuse gamma- 
ray background virtually excludes domains of antimatter in the visible Universe [if. 

Direct observations of luminous matter show that baryons constitute about a few percents 
of the total mass of the Universe. This gives a value of order 10 -10 for Ub/s, the ratio of the 
baryon number density to the entropy density. Precise measurements of the abundance of 
primordial light elements predicted in big-bang nucleosynthesis combined with the cosmic 
microwave background experiments restrict this ratio in the range ub/s ~ (5 — 10) x 10~ n 0. 

Many scenarios for baryon production in the early Universe have been proposed to explain 
this small baryon asymmetry, such as baryogenesis in grand unified theories, electroweak 
baryogenesis, leptogenesis, and Affleck-Dine (AD) baryogensis Affleck and Dine [H pro- 
posed a mechanism of baryogensis in supersymmetric models in which scalar quark and 
lepton fields obtain large vacuum expectation values along flat directions of the scalar po- 
tential. These coherent scalars or the condensate start to oscillate at a temperature of 
order 10 TeV when supersymmetry-breaking effects start to become important, and a net 
baryon number is developed and stored in the oscillating fields provided that C and CP 
symmetries are also violated. Subsequently, the scalar quark and lepton fields decay and 
produce a baryon asymmetry. In fact, the AD mechanism is too efficient and the resulting 
baryon asymmetry is usually too large. Several dilution processes have been considered to 
reduce the large baryon asymmetry to the observed value. They involve either introducing 
additional entropy releases after baryogenesis (e.g., by the decays of the inflaton [5J], massive 
field fluctuations 0, or the dilaton Q) or reducing the baryon production by invoking non- 
renormalizable terms In all of these studies, the scalar fermion fields have been treated 
as classical fields and obey the classical equations of motion. In this paper, we will be con- 
cerned with quantum fluctuations created from the dynamics of the Affleck-Dine scalar field 
and its back reaction effects. In particular, we study how these nonequilibrium processes 
affect the baryon production in the AD mechanism. 

It has been proposed that in a certain AD flat direction the condensate is not the state 
of lowest energy but fragments to form metastable or stable Q-balls, depending on the 
shape of the radiative correction to the flat direction For instance, this occurs along 
the flat direction with large mixtures of scalar top quark and light gaugino masses or in 
the if n L-direction. In this paper we estimate the baryon asymmetry assuming that the AD 
condensate does not lead to this type of Q-ball formation. Our presentation is organized as 
follows. The model of Affleck-Dine baryogenesis is described in Sec. II. We introduce the 
method of nonequilibrium field theory in Sec. III. It is then applied to the Affleck-Dine model 
in Sec. IV. We show in Sec. V how to incorporate nonequilibrium effects in the calculation 
of the baryon asymmetry, and present the results of our numerical calculation and discuss 
their implications in Sec. VI. Our conclusions are offered in Sec. VII. 



II. AFFLECK-DINE BARYOGENESIS 



Affleck and Dine |H have shown that, in a supersymmetric SU(5) grand unified model, 
there is a flat direction in the low-energy effective potential for the following set of vacuum 



states of the up (u), strange (s), and bottom (b) squark fields as well as the scalar muon 
field {p)\ 



(«§) = a, (ui) = v; (? 2 ) = a, (jl) = v; (ftj) = e*y/\v\> + \a\\ (1) 

with all other fields having vanishing vacuum expectation values. Here the superscript 
c denotes charge conjugation, the subscripts denote color indices, a and v are arbitrary 
complex numbers, and £ is real. A large baryon asymmetry can be generated by the decays of 
the condensate associated with these vacuum expectation values via baryon-number violating 
dimension-four scalar couplings which appear after supersymmetry has been broken. To 
illustrate how the mechanism works, Affleck and Dine considered a toy model with a single 
complex scalar field $ described by the Lagrangian density 

C = (d^) («) - m|$t$ - iX (V - $ t4 ) , (2) 

which contains the mass term with m% = Mj and the baryon-number violating coupling, 
A ~ eMg/Mg, where Ms is the effective supersymmetry breaking scale, e is a real parameter 
characterizing CP violation, and Mq is the grand unification scale. For small $, the theory 
has an approximately conserved current 

j, = % {fdfi - (0 M $t)$) ( 3 ) 

due to the approximate global U(l) symmetry: $ — > e j6l $. The corresponding charge will 
be referred to as the baryon number (this will be the case if, for example, $ represents a 
scalar quark field). 

Classically, in an expanding universe, $ obeys the equation of motion 

£ + aff f +m .. = «A.f, (4) 
where H is the Hubble parameter. With the initial conditions at t = t : 

$| f=fo =i$o and $| t=to = 0, (5) 

where $o is real and $ = d$>/ dt, it was found that the baryon number per particle at large 
times (t ^> mj 1 ) in either a matter-dominated or a radiation-dominated universe is given 
by 

r ~ A$g/m|. (6) 

Assuming m| = Mf and A = sM^/Mq, we have r ~ c^q/M^., which can easily provide 
a large initial ti-q/s and thus dilution processes have to be introduced to reduce it to the 
observed value. For example, taking e = 10 -3 , Ms = 10~ 16 Mp, Mq = 10~ 1 Mp, and 
$^ = 10 _3 Mp, where M P is the Planck mass, we find A = 10 33 and r ~ 10 4 . Note that 
the decay width of the condensate can be estimated as T$ ~ (a s /7r) 2 m|/|$| 2 [4], which is 
typically much smaller than the frequency of the oscillating scalar fields, i.e., T$ m$. 
Therefore, a net baryon number is developed and gets saturated in the oscillating scalar 
quark and lepton fields before they decay and produce a baryon asymmetry. 

The above consideration did not take into account the fact that a dynamical background 
field could generate quantum fluctuations which could in turn influence its evolution [§, 



loL 11, Substantial quantum fluctuations from the vacuum state may lead to a value 



of the generated baryon number density quite different from that obtained from purely 
classical arguments. In the following sections, we will study how the quantum fluctuation 
of the dynamical $ field and the back reaction on its evolution will affect the ratio r for 
the model of Eq. (J2J). Nonperturbative Hartree factorizations will be implemented along 
with the method of nonequilibrium quantum field theory to self-consistently take account of 
quantum fluctuations and back-reaction effects on the background field. To be completely 
general, we will not fix the values of m$, A, and <3>o, leaving them as free parameters. 



III. NONEQUILIBRIUM QUANTUM FIELD THEORY 

Before we discuss the Affleck-Dine baryogenesis, let us review the field theoretical methods 
for treating the quantum fluctuation and back reaction of a dynamical quantum field. The 
basic methods for studying nonequilibrium phenomena were developed many years ago by 
Schwinger and Keldysh [13[. In nonequilibrium situations we are interested in obtaining 
the equations of motion for the expectation values and correlation functions of the quantum 
fields in an evolving quantum state or density matrix. This is accomplished by tracking the 
time evolution of the density matrix, which corresponds to an initial valued problem for the 
Liouville equation. 

In the Schrodinger picture the evolution of the density matrix p is determined by 

p(t) = U(t,t )p(t )U- 1 (t,t ), (7) 
where U (t, to) is the time evolution operator: 



U(t, to) = Texp 



-i f dtH{t')\ (8) 



for a time-dependent Hamiltonian H. The symbol T means to take the time-ordered prod- 
uct. The expectation value of an operator O in the Schrodinger picture is given by 

/nVA _ Tr [Ujt^p^U^jt^O] 

{){) ~ Tr[p(t )] ' {9} 

and thus can be determined once the initial density matrix p(to) is specified. 

Consider the case in which the initial density matrix describes a system in equilibrium. 
When a perturbation is switched on at time to, the resulting time-dependent Hamiltonian 
can drive the initial state out of equilibrium. Just as in the usual perturbation theory, we 
write H = H + H int , where H is the unperturbed quadratic Hamiltonian and H int stands 
for the perturbations. For technical convenience of treating this initial valued problem, we 
can model the dynamics by the following Hamiltonian: 

H(t) = Q(to - t)H (to) + Q(t - to)[H (t) + H int (t)], (10) 

where 0(t — t ) is the Heaviside step function. For t < t , the unperturbed Hamiltonian 
H (to) is fixed at time to from which we can specify the initial thermal state of the system, 
and for t > t , the perturbation is switched on in H int . Before the perturbation is switched 
on, the system is assumed to be in equilibrium at a temperature T = 1/(3 and is described 
by the initial density matrix 

p{t ) = exp[-/3iJ (to)] ■ (11) 



The zero temperature limit for an initial vacuum state can be obtained by taking T — > 0, as 
we shall consider later. 

Notice that the initial density matrix can be expressed in terms of the time evolution 
operator along imaginary time in the distant past: p(t ) = U(—oo — i/3, — oo). This allows 
us to write the expectation value in Eq.fjUJ) as 

Tr[^(-oo-z/3,-oo)] ' [U) 

where we have inserted U(t , — oo) {7 _1 (£ , — oo) = 1 to the right of p(to) in the numerator, 
commuted U(to, — oo) with p(t ), and used the following composition property of the evo- 
lution operator: U~ x {t^ -co) ^"HMo) = [U(t,t ) U(t , -oo)] -1 = t/" -1 ^, -oo). Another 
insertion of 1 = £/ -1 (oo, t) [/(oo, t) yields 

= Tr [[/(-oo - jg, -oo) [/^(oo, -oo) [/(oo, t) O U(t, -oo)] 
Tr [U(— oo — z/3, — oo)] 

Although not apparent from the above expression, the time-dependent expectation value for 
O depends implicitly on £ through the perturbation H int which is switched on at time £ - 
Because there are both forward and backward [given by f/ _1 (oo, — oo)] time evolutions, and 
because the initial density matrix is written as the time evolution operator along imaginary 
time, one is led to consider the following generating functional 

Z\j+,r\ = Tr [[/(-oo - i(3, -oo) [/(-oo, oo, j~) [/(oo, -oo, j+)] (14) 

for deriving the nonequilibrium equations of motion for the expectation values and correla- 
tion functions of the quantum fields. 

In anticipation of the subsequent application to the Affleck-Dine model (J2J) for baryoge- 
nesis, let us consider the case of a scalar field p as an example. In this case the generating 
functional above has the following path-integral representation 

Z[j + ,j-] = J dipdtpxdipzj Vcp+Vcp'Vcp 13 exp (i J°° dt J rf 3 x {c [<^ + ] + j + (p + ]) 
xexp^-z J dt J d 3 x \C p>~ + j~yr]^expN J ^ dt J d 3 x |£ P P 

(15) 

with the boundary conditions: <^ + (x, — oo) = ^(x, — oo — i(3) = y?(x), <£> + (x, 00 ) = 
y?~(x, oo) = y?i(x), an d f"^, — oo) = ^(x, — oo) = <^2(x). The path integrals represented 
by the p + and p~ fields correspond to the forward and backward time evolution, respectively, 
while the one described by the p" field corresponds to the time evolution along the path 
parallel to the imaginary time axis. The boundary conditions are in fact a result of the trace 
as well as the bosonic nature of the operators. Notice that the path integral corresponding 
to evolution in imaginary time involves only the unperturbed quadratic Lagrangian £ [ see 
Eq. (jTUj) ]. while the other path integrals involve the full Lagrangian £ = C + £ int , where 
£i n t stands for the interaction Lagrangian that will be treated as a perturbation. 

The introduction of sources in Eq. (JT5)l allows one to obtain expectation values involving 
the quantum field operators by taking functional derivatives with respect to the sources as 
follows: 

P + - -i± , <p- - , (16) 



and then setting the sources to zero in the end. In fact, various Green functions can be 
obtained by taking functional derivatives of the generating function with respect to the 
appropriate sources. For instance, functional differentiation with respect to the sources j + 
and j~ can generate the time-ordered and anti-time-ordered Green functions, respectively. 

As for computing the real-time correlation functions of interest, there is clearly no need to 
introduce any source along the path parallel to the imaginary time axis since the Lagrangian 
along this path involves only its unperturbed part C . In the Schrodinger picture, the 
ensemble average we implement corresponds to an initial density matrix describing a state 
of thermal equilibrium with respect to the unperturbed Hamiltonian. Path integrals over 
Cq can be carried out exactly and we can obtain the generating functional in terms of the 
nonequilibrium propagators. The temperature dependence enters through the boundary 
conditions on the nonequilibrium propagators as we shall see later. Thus, the relevant 
generating functional for computing real-time correlation functions can be obtained as 



x exp < — 




C 



hit 



_6_ 



x) G ++ {x, x') ] + {x') - ] + {x) G + -{x, x') j-(x') 

r(x)G- + (x,x') 3 + (x')+r(x)G-(x,x , )r(x')}} , (17) 



where the nonequilibrium propagators corresponding to the forward and backward time 
branches are given by 



G ++ (x,x';M') 
G (x,x'; t,t') 
G + -(x,x';t,t') 
G _+ (x,x';t,0 



<^(x,tV+(x',t')> 

G > (x, x'; t, t') Q(t - t') + G<(x, x'; t, t') 6(t' 
^-(x,t)y9-(x',i')> 

G > (x, x'; t, e(t' - t) + G<(x, x'; t, t') Q(t - t' 

(^ + (x,t)^(x',t')) 
G < (x,x';t,t / ) = (^(x / ,t , )^(x,t)), 

(^-(x,t)^(x',t')) 



t) 



x, x 



t,t') = (tp(x,t)<p(x!,1?)). 



(18) 



Since time translational invariance is lost for a system that is out of equilibrium, the Green 
functions defined above will in general depend on the time t and t' separately. 

The functions G > and G < can be obtained from computing the 2-point correlation func- 
tions of the quantum field <p, which can be expressed in terms of the usual creation and 
annihilation operators with mode functions that are solutions to the homogeneous equa- 
tions of motion for the quadratic Lagrangian Cq. The explicit derivation for them will be 
shown in the next section. As a result of the assumed initial thermal state mentioned above, 
they obey the boundary condition 



G< (x, x'; -00, If) = G> (x, x'; -00 - i(3, if) , 



(19) 



which corresponds to the Kubo-Martin-Schwinger (KMS) condition at a large negative time 
before the pertubation is switched on |l5j |. 

In summary, the generating functional (|17|) can be used to obtain the Feynman rules that 
define the perturbative expansion for weak couplings. There are four sets of nonequilibrium 



propagators in the forms given in Eq. (|18j). Two sets of interaction vertices denned by the 
interaction Lagrangian £ int involve fields in the forward branch and those in the backward 
branch. Notice that there exists a relative minus sign difference between these two types of 
vertices. The calculation of the combinatorial factors is the same as in the usual quantum 
field theory. This constitutes the basic tool for our studies. 



IV. COMPLEX SCALAR FIELD IN AN EXPANDING SPACE-TIME 

Since the interaction rate for scattering processes between the scalar fermion fields and 
the thermal bath in the early Universe is found j|4] to be of the same order of magnitude 
as their decay rate Tg 1 (see Sec. II), the scalar fermions are essentially decoupled from the 
bath for a time scale much longer than the oscillation period of the scalar fields: mj 1 <C 
t <C Tg 1 . We can therefore apply the method of nonequilibrium quantum field theory with 
zero temperature to study Affleck-Dine baryogenesis by considering the model of Eq. (J2J). 

The relevant action can be written as 

S = Jd 4 x [ - m%&$ - %\ ($ 4 - $ t4 ) ] . (20) 

The background geometry is governed by the spatially flat Robertson- Walker metric, 

d 2 s = d 2 t-a 2 {t)d 2 x, (21) 

where a(t) is a scale factor. It will be convenient to introduce the conformal time rj defined 
as drj = a~ 1 dt and the conformally rescaled field, \ = a$, whereby the above action can 
be written as that of a complex scalar field in flat space-time with a time-dependent mass 
term. To see this, we express the complex scalar field \ in terms of its real and imaginary 
parts: % = (Vv^) (xi + %)• The action becomes 



S = J dr)d 3 x£ [xi,x 2 ] 



drj d x 



Vi 2 - \ (vXi) 2 + \x2 - \ (VX2) 2 -U(xi,X2) 



(22) 



2 2V V/V V 2 
where \\ = dxi/dr] and the new potential 

U( X i, X2) = \m 2 x {r,) (xl + xl) ~ 2A X iX2 (xl ~ xl) (23) 

has the time-dependent mass term with m^rj) = m\a 2 — (a" /a). 

Since we are interested in the effects of the quantum fluctuations of the scalar fields xi 
and Xii it wm b e useful to express each of them in terms of its expectation value with respect 
to an initial density matrix that will be specified later and the fluctuation field around this 
average value: 

xf{*,v) = xKv) +xt(*,v) , (xt(*,v)) =xi(v) ; 

xf(*,v) = xl(v) + xf(*,v) , (xtfav)) = Xiiv) ■ ( 24 ) 

We have listed the fields belonging to the forward (+) and backward (— ) time branches 
explicitly. The reason for shifting the (±) fields by the same configuration is that the field 



expectation values enter in the time evolution operator as c-number background fields, and 
the evolution forward and backward are considered in this background [9j. Note that Xi 
and X2 are time-dependent for the nonequilibrium problem. We will implement later the 
corresponding tadpole conditions, {xf (x, 77)) = and {xf (x, 77) ) = 0, to derive the equations 
of motion for these field expectation values. Expanding the action in terms of the fluctuation 
fields xi(x, 77) and X2(x, 77), the Lagrangian density becomes 



C 



xl + xt,xl + xt 



c 



xixl 



+ c 



xt,xt 



- ^|x?"+^)Xi-2A[3(x?) xl~{xl 

- xt {xl"+m 2 x (v)xl + 2X [3x?(^) 2 -(x? 

+ ^xlxl(xt 2 -xi 2 ) + ^xtxt 

+ 2 A (xlxt 3 - xlxV + ZxlxVxt - -WM\t 2 



xtf-(xf 



(25) 



and similarly for £ [xi + Xi , xl + X2 

We now introduce a Hartree factorization to take account of the quantum fluctuations. 
The Hartree factorization can be employed for both the (+) and (— ) fields as follows: 



xi ~ 

X 2 iX2 

X\X2 



HxDxi 



xl ~ 
xixl 



Hxl)x2, 



2<XiX2>Xi + (%i)X2, 
3(XiX2)x? + 3(x?)XiX2, Xixl ~ 



2(xiX2)X2 + (x 2 >Xi 5 
3(XiX2)X2 + 3(x1)XiX2 



(26) 



where the one-loop corrections to the 2-point Green functions of xi an d X2 are cancelled by 
the counterterms introduced above. We would like to point out that the Hartree approx- 
imation we adopt here is uncontrolled in the model of Eq. (J2J). It will be controlled, for 
example, in alternative models in which the scalar fields are in the vector representation of 
the symmetry group O(N) as the Hartree approximation is then equivalent to the large- N 
approximation. Our justification for using this approximation here is simply based on the 
virtue that it provides a nonperturbative framework to treat the quantum fluctuations by 
solving the evolution equations self-consistently. 

After applying the above factorizations, the nonequilibrium Lagrangian density can be 
written as 



where 



C 

Co 



Xi + xt,x°2 + xt 



c 



Xi + Xi , xl + X2 



Xi,X2 



Xl,X2 



Co 



Xi ,X 2 



c 



ltlt 



Xl",X2 



c 



int 



Xl ,X2 



(27) 



2*i 



1 

2 (yxi 



+ 



\xi' 2 - \ (vx^) 2 - lM 2 X2 (v)xi 2 , 



(28) 



with 



nr 



rn 



\{ri) - 12Ax?X° - 12A(xiX2> , 
:(r7) + 12A X ?X2 + 12A(xiX 2 ), 



(29) 



and 



Ant [xi, xt\ = -an{n)xi - a 2(v)xt + ®i2(jti)xtxi , (30) 

with the couplings 

ai(v) = X?" + mlWxl - 2A [3 (xtf xl ~ (xtf] ~ ^X° 2 (<X?> - (xl>) ~ 12Ax?<XiX2> , 
a 2 fo) = xl" + m£fa)x§ + 2A [s X ? (x^ - (xi)*] - 6A X ? (<XI> - (x1>) + 12A^<XiX 2 > , 



atniv) = 6A 



Xi 



xl 



+ 6A ((xl) - (x1>) ■ 



(31) 



We have used the fact that correlation functions of fields evaluated at the same space-time 
point in the forward and backward time branches are equal. This can be seen from Eq. (I18|) 
by identifying (x, t) with (x',t'). In addition, due to the spatial translational invariance, 
the correlation functions depend only on time and should be understood as: (xf 2 (x,r])) = 

(Xi 2 (*,v)) = (xl)(v) (i = 1,2), (x^(x,r])x 2 f (x,r/)) = (xT(x, v)x.2 (x, rf)) = (xiX 2 )(^)- 

In Eq. (j27|) . £o is the unperturbed quadratic Lagrangian with respect to which the 
nonequilibrium Green functions will be defined and the interaction Lag rangian £; n t con- 
tains the linear terms in \i as we U as the X1X2 terms which will be treated as perturbations. 
The presence of the nondiagonal X1X2 terms implies that these fields do not form a good 
basis to construct the nonequilibrium propagators in a perturbative expansion. One may 
therefore consider invoking the Bogoliubov transformation which is a canonical transforma- 
tion in the sense that it leaves the measure of the path integral invariant. In this case, the 
transformation is a two-dimensional field rotation between the old and new basis fields with 
one parameter specified by the rotation angle. The angle for the field rotation can be deter- 
mined so as to rid the nondiagonal term in the potential energy of the Lagrangian. In the 
presence of the time-dependent background fields, it turns out that the Bogoliubov trans- 
formation acquires a time-dependent rotation angle, which in turn gives rise to a new kind 
of unwanted nondiagonal terms in the kinetic energy. Thus, it seems that the Bogoliubov 
transformation fails to diagonalize the Lagrangian in the presence of the time-dependent 
background fields. To proceed, we can further approximate the term X1X2 by its mean value 
(XiX2)- In other words, we write 



X1X2 = (X1X2) + (X1X2 - (X1X2)) 



(32) 



and treat the terms in the parenthesis as a perturbation. 

The fluctuation fields can be decomposed in terms of their Fourier modes with time- 
dependent mode functions as 



Xi&v) 



ikx 



y= Z [ fl i>k fiM + 4,-k /i%kfa) 



ikx 



1,2. 



(33) 



where a,^ and a\ _ k are, respectively, the annihilation and creation operators for the fluctu- 
ation field Xi with momentum k. The equations of motion for the mode functions /j^ can 
be read off from the quadratic Lagrangian £ i n Eq. (|25|) as 



A 2 

di] 



2 +k 2 +M 2 Xi ( V ) 



fiM = 0, i = l,2. 



(34) 



where k 2 = k ■ k. The nonequilibrium propagators can also be expressed in terms of the 
mode functions: 



G t*(v,v 

Gf k (v, V 



GUv, rf) 6(77 - i) + G< k (7?, rf) Q{rf - 77) 
G> k (v, rf) 8(77' - 77) + G< k (r/, rf) 6(r/ - rf) 
GUvW), 
GUV,V'), 

|d 3 x(x,(x,r ? )x l (0,r ? '))e- kx 



(35) 



where the correlations of the fluctuation fields are taken with respect to a vacuum state that 
is annihilated by the cij k of Eq. (l3~3*|) . 

The equations of motion for the expectation values of the fields can be obtained from 
the tadpole conditions. From the condition (xi~(x, 77)) = 0, the lowest-order contributions 
to the equation of motion are obtained from the terms linear in xt m the interaction La- 
grangian ()3(Jj). We find 



drj -ui{r))(xtxt) + aci{r}){xtxi) = 0, 



(36) 



where the spatial arguments have been suppressed. Since the correlation functions (xtxt) 
and (xiXi ) are linearly independent, the above equation leads to the equation of motion 
for Xv 

XT + mjfo)*? ~ 2A [3 ( X ?) 2 X° 2 ~ (xl) 3 ] ~ ^X% ((xl) - (xl)) ~ ^xliXiX*} = , (37) 

where Eq. (13 1JI has been used. It is straightforward to see that implementing the condition 
(Xi"( x ? v)) = w i n l ea d to the same equation of motion for Xv Following the similar 
procedure with the conditions (% 2 (x, 77)) = 0, we find the equation of motion for x\ to be 



xX +m 2 ( V )x° 2 + 2\ 



3X? (Xl 



xl 



6Ax? {(xl) - (Xl)) + 12A X °(xiX2> = . (38) 



The quantum fluctuations (x\) and (%|) can be determined self-consistently and have the 
form 

<X?>fo) = / 7^)3 I I 2 , (39) 

where the infinite volume limit has been taken. On the other hand, the term (X1X2) vanishes 
at tree level and one has to include an interaction vertex from the perturbations to obtain its 
loop corrections. To lowest order in perturbation theory, the vertex QiX (xl) 2 ~ (x°) 2 X1X2 
is included and one finds that 



(XiX2)(v) 



d 3 k 

(2tt) 3 



(xi,Mx2,-M) 



6t\ drj \ XiW) 



6i\l dri\ xi(v') 



-(xiM(v)xi,-k{v'))(xt-Mx2,k(v')) 



(40) 



As mentioned earlier, the 7] dependence arises from the perturbations being switched on 
at time r) . We have shown above only the calculation (second equality) for the X1X2 
correlation. However, as we have stated after Eq. (}3T|) . the final result is also applicable to 
(xTxT)- m terms of the mode functions, this correlation function becomes 



(XiX2)(v) = 6iA / drf { XiW) 



xl(v n 



fi,k(v)fi,k(v')h,k(v)f2,k(v') - c - c ] ■ (41) 



Equations (jMJ), (jSZj), (jHHJ) and (jUJ) form the full set of coupled equations from which 

we will solve for the field expectation values and the mode functions. 



V. NONEQUILIBRIUM EFFECTS AND BARYON ASYMMETRY 



The proper normalization to measure baryon asymmetry is given by the particle number 
density. Let us consider the initial state at time to to be specified by the adiabatic modes 
for uncoupled harmonic oscillations defined in cosmic time under Hartree approximations. 
The corresponding initial particle number density can be expressed as 



n(t ) 



d 3 k ( a 3 (t c 



(2tt) 3 {2W 1M (t 
a 3 (t ) 



$i,k(*o)$i,-k(*o) + Wjyto) $i,k(to)$i,-k(*c 



2W 2 , k (t ) 



$2,k(*o)$2,-k(*o) + WL(t ) $2,k(t )$2,-k(tc 



(42) 



where $ = ($i + i$2)/v / 2- The instantaneous frequencies for the real and imaginary part 
of the complex scalar field $ can be gleaned from Eqs. (J3lj) and (}29|) to be 



W 2)k (t ) 



k 

— + m%- 12Aa- 2 (t ) (xi(*o)x!i(to) + (xiX2>(*o) 
k 2 

^— - +m% + 12Aa- 2 (t ) (xi(*o)xjj(*o) + (xiX2>(to) 



(43) 



The expectation value of the number operator with respect to an initial vacuum state as 
defined after Eq. evolves in time and has the following form in the Heisenberg picture: 



(n)(t) 



d 3 k \ a 3 (t ) 



(2tt)3 [2W 1)k (t 
a 3 (to) 



($i, k (t)$i,- k (t)> + WlJto) ($i, k (t)$i,_ k (*)) 



2W 2>k (t c 



($2,k(t)$2,-k(t)> + m, k (t ) ($2,k(t)$2,-k(*)> 



(44) 



In terms of the conformally rescaled fields, the expectation value of the particle number 
density can be written as the mean value constructed from the expectation values of the 
quantum fields plus its quantum fluctuations: 



(h)(r]) = n(r}) = n (r]) + n q (ry) 



710(77) 



n q (rj) 



a 3 M 



2Wi, (T/o)a 4 (77) 



A. 

dq 



xl(v) 



a!{rj) d 
a{rj) dr] 



+^i,k(r?) \x°i(r)) 



+ (xi X2 ; tti,k(v) -> fi 2,k(^)) 



d 3 k 



a 3 (770) 



(2tt) 3 i2Wi, k (7/ )a 4 (77) 
+^i,k(r/) (Xi,k(r/)xi,_k(??))] 



(xUWxi.-kW) "44 f#(Xi,k(r/)xi,-k(r/)) 



0(77) \<i?7 
+ (Xi X2 ; ^i,k(r/) -> ^ 2 ,k(??)) 



where 



^,k(r/) = Wi(7 /0 )a 2 (r / ) + 



a 2 (77) 



(45) 



(46) 



and 77 = 77(^0) is the initial conformal time. 

The correlation functions can be expressed in terms of the mode functions as in Eq. (|39j): 



(47) 



In the weak coupling limit, and with appropriate choices for the initial expectation values 
of the quantum fields, we may approximate 



Wi, k (t ) « W 2 , k (t ) « W k (to) ; W k (t c 



A' 2 



a 2 (t ) 



+ rn% 



(48) 



This approximation simplifies the expression for the baryon number density which is given 
by the current in Eq. (JHJ) : tt-b = — $^$). In parallel with the discussion above for the 

particle number density, we find 



n B (t) 



a 3 (t ) 



d 3 k 

(X) 3 



$l k (t)$ k (t)+6_ k (t)4(t)-$L k (t)$ k (t)-$_ k (t)$ k (t) . (49) 



After converting to conformal time, the expectation value of the baryon number density is 
found to be 



n B (v) = {nB,)(v) 



a3M f Mv)x?(v)-x°i(v)xl'(v) 



a 3 (r]) 
a 3 (770 



d 3 k 



a 3 (r]) J (2vr) 3 



(xi,k(?7)X2-k(^)) - (Xi,k(?7)X2-k(^)} • (50) 



Like the particle number density, it can be separated into a mean value plus its quantum 
fluctuations. Just as in the earlier computation for the term (X1X2); the quantum fluctua- 
tions can be computed by introducing a vertex 6iA [(x?) 2 — (X2) 2 ] X1X2 at the lowest order, 



which yields 

a 3 (r] 



n B (r}) 



a 3 (rj) 



xl(v)x°i (v) -Xi(v)x°2 (v) 



+6ia% Vo )X / drj' { X °iW) 
'no 



c.c. 



xlin')] }{[fUv)ftM)f2Av)fiM 
fiAv)ft*(v')fUv)fi*(v') - cc] } . (51) 



The baryon number density can be evaluated once the mode functions and field expec- 
tation values are determined by the method described in Section IV. In order to ensure the 
canonical quantization rules, the mode functions must satisfy the relation 



(52) 



which follows from the mode equations (|34p. With the initial states described by the adi- 
abatic modes for uncoupled harmonic oscillations, the mode functions have the following 
initial values at the conformal time f] : 



fi,*(Vo) 



a(r}o 



2a 3 MW k (7? ) 



fU(Vo) 



2a 3 (%)W k ( % ) 



— i 



a 2 M WkM 
2a 3 (r /0 )W k (r /0 ) 



(53) 



As expected, these initial conditions guarantee that the baryon and particle number densities 
have zero quantum fluctuations initially. 

Consistent with the initial conditions (J5)l. the initial expectation values of the conformally 
rescaled fields are 

X%Vo) = , Xi(Vo) = ; xliVo) = V2a( Vo ) $ , X °2 (Vo) = V2a'( Vo ) $ . (54) 

Together with (|53J) . these initial conditions guarantee that the initial baryon asymmetry is 
zero: n-Qirjo) = 0. 

In the next section, we will present the results of our calculation that assumes a radiation- 
dominated epoch in which the scale factor follows a(i]) = m$r]. We find that the results are 
similar in a matter-dominated universe and will not present explicit results for this case. 



VI. RESULTS 



Before we describe the numerical solutions to the coupled equations, let us try to under- 
stand the issue of ultraviolet divergence associated with the loop integrals in these equations. 
The large-fc behavior of the mode functions can be obtained from WKB type solutions to 
the mode equations: 



AM 



a(r) ) 



2J2a3(r ?0 )W k (r /0 ) 



;i + T)^k(^) + a - 7)A,kfo) 



(55) 



where 



(56) 



and 

n„M = -* - ±MUn) - ~MIM + o (1) + ■ ■ ■ . (57) 

The parameter 7 that appears in Eq. (|55j) can be determined from the initial conditions (|53j) 
on the mode functions. It has the large-fc expansion 

7 = -l-£(^WiV-. <«) 



which secures 

«^-^H>»-3^) + °(f) + '''- < 59 » 

These large-fc behaviors lead to quadratic and logarithmic divergence for the loop integrals 
as follows: 

itm = ^ - A i M lM - ^rr) ln - +finite ' (60) 

where A is the ultraviolet momentum cutoff and k is the renormalization scale. Let us 
recall that the Affleck-Dine model under consideration is used as an effective field theory to 
describe the relevant low-energy physics from a supersymmetric grand unified theory. A is 
therefore a physical cutoff which can be set at the scale of grand unification, A m Mq, while 
k can be chosen to be the supersymmetry breaking scale, k ~ Mg ~ Since the relevant 
terms involved in the equations for the field expectation values are of the form ((Xi) ~ 
only the logarithmic dependence on the cutoff A remains. For weak couplings, such a mild 
cutoff dependence does not present any difficulties for our numerical study. Using Eqs. (jS3J), 
and (|57jl. the quantum fluctuations (X1X2) prove to be ultraviolet safe. 

Let us turn now to the numerical results of our study. Figure 1 shows the evolution of 
the expectation values of the quantum fields Xi(v) an d Xi.i 7 !) m units of m$/vA, the scale of 
which is set by the dimensionless parameter O = v^A$o/ m * = 1CT 2 . According to Eqs. (|S7j) 
and (|3H|). the evolution of the field expectation values are determined by the fluctuations 
((xl) - (X 2 2 ))(v) and (xiX.2)(v), which are displayed in Fig. 2. 

For an extremely weak coupling, e.g., A = 10~ 20 , Fig. 2(a) shows that the quantum 
fluctuations are very small and Xiiv) an d X^i 1 !) follow essentially their classical trajectories 
which correspond to oscillations with time-dependent frequencies. For a larger, but still 
weak, coupling (e.g., A = 1CT 3 ), the fluctuations become more significant (see Fig. 2(b)) and 
the amplitudes of the oscillations of Xi an d X2 become damped, as can be seen from Fig. 1. 
From Eq. (JHTj) and the initial conditions (p>fj) . we can estimate the order of magnitude of 
Xi(v) to be 

xS(l)«Aif, (61) 

which is consistent with Fig. 1. 

The time evolution of the baryon asymmetry, as measured by the ratio of the baryon 
number density to the particle number density, r(rj) = n^{rj) j 'n{j]) , is shown in Fig. 3. As 
stated earlier, the initial baryon asymmetry is zero due to the initial conditions in Eqs. (J53|) 
and (JSH). For very weak couplings, the evolution of r is determined by the classical dynamics 





FIG. 2: Time evolution of (xf ) — (smooth curves) and (X1X2) (oscillating curves) for 4>o = 10 2 
and (a) A = 1(T 20 and (b) A = 10~ 3 . Note that we have used the rescaled fields Xi = v~\Xi/m$. 




FIG. 3: Time evolution of the ratio of the baryon to particle number density, r(rj) = n-B(rj) / n{rj) , 
for (p = 1(T 2 and (a) A = 1(T 20 and (b) A = 1(T 3 . 




FIG. 4: (a) Time evolution of the ratio of the total particle number density to its classical value 
for cj)Q = 10~ 2 and A = 10~ 3 . (b) Same as (a) but for the baryon number density. 
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FIG. 5: Ratio of the baryon to particle number density at large times for various initial values of 
$o and couplings A. 



and is depicted in Fig. 3(a). It starts from zero and decreases toward an equilibrium value 
r(oo) ~ — 1000- This saturation value can be understood as follows: 



r oo 



n B (oo 



n oo 



{xlf 



{xl 



ml 



(62) 



where the approximation for Xi in Eq. (JrTTj) has been applied to compute the baryon number 
density ne, but neglected in the number density n as compared with that of x%- The negative 
baryon asymmetry should not be a concern as the sign can be fixed by the baryon number 
assignment for the field The alert reader will recognize that the estimate of r(oo) above 
is nothing but the result quoted in Eq. (JHJ). 

As can be seen from Fig. 3(b), an interesting phenomenon occurs for larger couplings when 
the quantum fluctuations become large, typically of the order of the field expectation values. 
The growth of quantum fluctuations is a consequence of parametric amplification [l2l . IT^ 
induced by the oscillating mean fields. In this case the part of the ratio r{rf) generated from 
the quantum fluctuations turns out to dilute the mean-field contribution, and therefore drives 
the ratio toward a much smaller equilibrium value: |r(oo)| pa lO^A^/ml,. It thus appears 
that, under appropriate conditions, the effects of the nonequilibrium quantum fluctuations 
can wash out the baryon asymmetry generated from the classical dynamics. The reason for 
this may be understood as follows. The set of coupled equations, Eqs. (J33|) . (pTTjl . (jSBJ), (f3TJj) 
and (|41jl. which we solve to obtain the baryon asymmetry, involves fluctuation terms with 
nonzero baryon number such as ((x?) — (xi)) an d (XiXz)- such, the amount of baryon 
asymmetry that can be generated is limited by the Hartree approximation implemented 
in these equations. On the other hand, as a result of the nature of the couplings under 



consideration, it is possible to create a large particle number density from the fluctuation 
term ((xi) + which has nonzero particle number but zero baryon number under the 

global U(l) symmetry. This can be seen in Fig. 4, which shows that for a sufficiently large 
coupling the total particle number density deviates substantially from its classical value, 
while the baryon number density remains close to its classical value. A similar plot for a 
much weaker coupling, e.g., A = 1CT 20 , would show that both the particle number density and 
the baryon number density remain essentially at their respective classical values. We are thus 
led to the conclusion that the ratio r can be reduced significantly through nonequilibrium 
quantum fluctuations with a proper choice of the couplings. 

Figure 5 shows the saturation values of r as a function of the coupling A for a wide range 
of values of $ - Plotted are equal O contours, where (f) = y/X^o/mq,. For a given contour, 
the saturation value for small values of A is given by the classical result, |r(oo)| ~ 1O0q, 
which corresponds to the horizontal part of the curve. When A becomes larger than about 
1O _6 0O) the nonequilibrium effects start to operate and significantly diminish the saturation 
ratio. For a wide range of parameters, this can result in r saturating toward the observed 
value of the baryon asymmetry, as shown by the right end-points of the contours. 

As an illustration, let us recall the sample values of the model parameters cited after 
Eq. ©, which corresponds to the left end of the 0o — 10~ 2 contour in Fig. 5. Now let 
us entertain the possibility that the supersymmetry breaking scale may be as high as the 
grand unification scale, i.e., Ms ^ Mq- We then find A % 1CT 3 and the baryon asymmetry 
is reduced by the nonequilibrium dynamics to the present observed value, |r(oo)| ~ 10~ 10 , 
corresponding to the right end of the 0o — 10~ 2 contour in Fig. 5. 

VII. CONCLUSIONS 

We have studied the effects of nonequilibrium dynamics in the Affleck-Dine mechanism of 
baryogensis. Within the context of the model (j2J), we find that, for very weak couplings, the 
ratio of the baryon to particle number density is approximately 10A$Q/m|, as dictated by the 
classical dynamics, but for larger couplings it is driven by the nonequilibrium effects toward a 
smaller value given by \r\ « 10~ 5 A$Q/m|. Interestingly, in some cases these nonequilibrium 
effects can reduce the baryon asymmetry generated in the classical Affleck-Dine mechansim 
to the present observed value without having recourse to additional dilution processes. This 
encouraging result suggests that nonequilibrium dynamics can play a significant role in 
the evolution of the early Universe. It is therefore important to explore other possible 
manifestations of these nonequilibrium effects. The methodology developed and presented 
in this paper will be useful for such studies. 
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